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Abstract. - We derive the proper form of Virial theorem for a system of rotating self-gravitating 
Brownian particles. We show that, in the two-dimensional case, it takes a very simple form that 
can be used to obtain general results about the dynamics of the system without being required 
to solve the Smoluchowski-Poisson system explicitly. We also develop the analogy between self- 
gravitating systems and two-dimensional point vortices and derive a Virial-like relation for the 
vortex system. 
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Introduction. — The theory of Brownian motion is 
a fundamental topic in physics since the pioneering works 
of Einstein and Smoluchowski [1]. Brownian theory has 
applications in many other area including chemistry, bi- 
ology and finance. The study of Brownian particles in 
interactions was first developed in connection to colloidal 
suspensions, inhomogeneous fluids, and supercooled liq- 
uids exhibiting glassy features [2-4]. In that case, the 
interaction is short-ranged. The study of Brownian parti- 
cles with long-range interactions is more recent. For ex- 
ample, the case of self-gravitating Brownian particles has 
been studied by Chavanis and Sire in a series of papers 
(see [5] and references therein). For this system, the mean 
field approximation is exact in a proper thermodynamic 
limit N — > +oo with r\ = (3GMm/ R d ~ 2 fixed (where d is 
the dimension of space) and it leads to the Smoluchowski- 
Poisson system. The Smoluchowski-Poisson (SP) system 
displays many formal analogies with the Keller-Segel (KS) 
model [6] of chemotaxis for bacterial populations in biol- 
ogy (see [7] for a description of this analogy). In recent 
years, the SP system and the KS model have been the 
object of intense research from physicists (see, e.g., [5] 
and references therein) and applied mathematicians (see, 
e.g., [8] and references therein). 

For self-gravitating systems, the Virial theorem plays a 
very important role [9]. For a two-dimensional (2D) self- 
gravitating Brownian gas, it has been shown in [5, 7, 10] 
that the Virial theorem takes a very simple closed form. 



In this Letter, we shall extend this relation to the case 
of a rotating self-gravitating Brownian gas. So far, there 
is no result on the rotating case. From the Virial theo- 
rem, we can deduce general properties of the rotating SP 
system without being required to solve these equations ex- 
plicitly. In this Letter, we focus on the basic results and 
make simplifying assumptions: (i) mean field approxima- 
tion; (ii) single species system; (iii) overdamped dynamics. 
The general case, relaxing these assumptions, will be con- 
sidered in a companion paper (in preparation). We shall 
also develop the analogy between self-gravitating systems 
and 2D point vortices initiated in [11] and derive a Virial- 
like relation for the vortex system. 

Virial theorem for a rotating self-gravitating 
Brownian gas. 

The Virial theorem. We consider a single species sys- 
tem of self-gravitating Brownian particles with individual 
mass m rotating with angular velocity £1. We work in the 
rotating frame. In the mean field approximation, and in 
the strong friction (overdamped) limit £ — > +oo, the dy- 
namical evolution of the density p(r, t) is described by the 
Smoluchowski-Poisson system [5,12]: 



dp 
dt 



V 



(Vp + P V$ e// ) 



A$ = S d Gp, 



(1) 



(2) 



Pierre- Henri Chavanis 



where 



* e //(r,i) = 



(ft x if 



(3) 



is the effective gravitational potential accounting for in- 
crtial forces and Sd is the surface of a unit sphere in d 
dimensions (in this paper, d = 3 or d = 2 when H ^ and 
d is arbitrary when ft = 0). For the sake of generality, 
we consider the generalized Smoluchowski-Poisson (GSP) 
system [12] with an arbitrary barotropic equation of state 
p = p(p) but, for ordinary self-gravitating Brownian par- 
ticles, the equation of state is the isothermal one 



(4) 



Since the temperature is fixed, we are describing a dissi- 
pative gas in the canonical ensemble. 

The Lyapunov functional associated with the isother- 
mal SP system is the free energy F = Ej — TS B where 
Ej = | J p$ dr — i J p(ft x r) 2 dr is the Jacobi energy and 
S B = —k B f — In — dr is the configurational Boltzmann 
entropy. For the GSP system, the Lyapunov functional 
is the generalized free energy F = Ej + J p J p ^fr- dp' dr 
[12]. Using Eqs. (l)-(2), we easily derive the canonical 
//-theorem 



At equilibrium, F = 0, implying 

Vp + pVQeff = 0, 



(5) 



(6) 



which is the condition of hydrostatic equilibrium in the 
rotating frame. For the isothermal equation of state (4), 
we obtain, after integration, the mean field Boltzmann 
distribution 



p(r) = J 4e-/ 9TO *«//« > 



(7) 



where (3 = l/(k B T). It can be shown that this distribution 
is dynamically stable with respect to the SP system if 
and only if (iff) it is a (local) minimum of free energy 
at fixed mass. In that case, it represents the statistical 
equilibrium state in the canonical ensemble [5]. In the 
microcanonical ensemble, the statistical equilibrium state 
is a (local) maximum of entropy at fixed mass, energy and 
angular momentum [13]. The critical points (canceling the 
first order variations of entropy) are also given by the mean 
field Boltzmann distribution (7). Therefore, all the results 
obtained in the following at equilibrium, that exclusively 
use the form of Eq. (7), are valid both in the canonical 
and in the microcanonical ensembles 1 . 



lr The inequivalence between the microcanonical and the canoni- 
cal ensembles appears only when we consider the thermodynamical 
stability of the system which is related to the sign of the second 
order variations of entropy or free energy [14,15]. 



The equilibrium scalar Virial theorem for self- 
gravitating Brownian particles can be derived as follows. 
Taking the scalar product of Eq. (6) with r, and using the 
identity §V((ft x r) 2 ) = -ft x (ft x r) = ft 2 r - (ft ■ r)ft 
we obtain 



Vp • r dr - W u - /ft 2 + 7 4fc ft fc ft 4 = 0, 



(8) 



where we have introduced the tensor of inertia 7y = 
/ pXiXj dr (its trace is the moment of inertia I = J pr 2 dr) 
and the Virial W u = - / pr ■ V$dr. For d ^ 2, W u = 
(d— 2)W where W = \ J p$ dr is the potential energy and 

for d — 2, Wu = —^jf- [5]. Integrating the first term by 
parts, we obtain 



d pdr + W u + /ft 2 - / lfe ftfeft J = dPV, 



(9) 



where we have introduced the average pressure of the sys- 
tem on the boundary of the domain P = § pr ■ dS 
(if the pressure is uniform on the boundary: p(r) = pt,, 
then P = pt). For an unlimited system, PV — > 0, pro- 
vided that the pressure decreases sufficiently rapidly with 
the distance. For the isothermal equation of state (4), we 
obtain 



dNk B T + W u + /ft 2 - / 4fc ft fc ft l = dPV. 



(10) 



The out-of-equilibrium Virial theorem is obtained by 
taking the time derivative of the moment of inertia and 
using the GSP system Eqs. (l)-(2). After straightforward 
algebra, we obtain 

^i = djpdr + W rl + /ft 2 - /, fe ft fc ft 4 - dPV. (11) 

For an (isothermal) Brownian gas 

^1 = dNk B T + W vl + /ft 2 - /.fcftfeft, - dPV. (12) 

The two-dimensional case. For a 2D system of self- 
gravitating Brownian particles, the Virial theorem takes a 
particularly simple form. In that case, ft = ftz, <& e ff = 
$ - ift 2 r 2 and V(ftr 2 ) = 2ftr. The scalar Virial theorem 
(11) becomes 



1 f GM 2 

-£I = 2 pdr — + /ft 2 - 2PV. 

For the isothermal equation of state (4), we obtain 

i CM 2 

-£/ = 2Nk B T + /ft 2 - 2PV. 

2^ 2 

Introducing the mean field critical temperature 

GNm 2 



knTr. = 



it can be rewritten 
1 



-H = 2Nk B (T - T c ) + m z - 2PV. 



(13) 



(14) 



;i5) 



(16) 
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This relation is the main result of this Letter. It extends with 
the Virial theorem given in [5,7,10] for non-rotating sys- 
tems (O = 0). At equilibrium, 



PV = Nk B (T-T c ) + -m 2 



This can be viewed as the global equation of state for a 2D 
rotating self-gravitating isothermal gas. As explained pre- 
viously, it is valid both in the canonical and in the micro- 
canonical ensembles. For a non-rotating system (f2 = 0), 
it reduces to 



PV = Nk B (T-T c ). 



(18) 



This equation of state has been previously derived by 
Salzbcrg [16], Katz & Lyndcn-Bcll [17] and Chavanis 
[7,10], using different methods. Since P > 0, this rela- 
tion shows that equilibrium states can possibly exist only 
for T >T C . In an infinite domain (P = 0), they can possi- 
bly exist only at the critical temperature T = T C (see [5,7] 
for more details and explicit solutions). For a rotating sys- 
tem in an infinite domain (P = 0), the equilibrium Virial 
theorem reduces to 



i/fi 2 = Nk B (T c -T). 



(19) 



Since in 2 > 0, this relation shows that equilibrium states 
can possibly exist only for T <T C . 

In an infinite domain (P = 0), the Virial theorem (16) 
for a rotating self-gravitating Brownian gas reduces to 



1 



-n 2 I = 2Nk B {T -T c ). 



(20) 



It is worth noticing that this equation is closed, which is 
not the case for the Virial theorem of Hamiltonian sys- 
tems [9]. This simplification for Brownian systems is due 
to the strong friction limit and to the two-dimensional as- 
sumption. Solving this equation, we find that the moment 
of inertia evolves like 



I(t) 



1(0) 



2Nki 



-(T c -T) 



2Nk B , 



(21) 



This relation shows that the system will experience col- 
lapse or evaporation depending on the sign of the quantity 
in bracket. Let us introduce the new critical temperature 



i(o)n 2 

2Nk B 



(22) 



This critical temperature depends on the initial value of 
the moment of inertia and on the angular velocity. It can 
be written 



r (fi) = T c (i--5 



(23) 



GM 



21(0) 



1/2 



(24) 



(17) 

We note that T < T c . In terms of the critical temperature 
(23), the evolution of the moment of inertia is given by 



I(t) = 



2Nk £ 
~ti 2 ~ 



(T c -T) + (T-T )e — 



For small times, we have 
ANk B , 

J(*)~I(0) + — ^*(T-T )t, 



(t-0). 



(25) 



(26) 



For T > To, the moment of inertia I(t) increases and 
tends to +oo as t —* +oo. Since I(t) = M(r 2 )(t), this 
corresponds to an evaporation process. For T < T$, the 
moment of inertia I(t) decreases and becomes zero in a 
finite time 



* - J_i 

tend - HI 



T c -T 

Ta-T 



(27) 



This corresponds to a finite time collapse. At t = t en d, the 
moment of inertia 1 = 0, implying that the system forms 
a Dirac peak containing the whole mass M. For T = To, 
the moment of inertia remains constant for all times: 



To) = 1(0). 



For a non rotating infinite system (Q = P 
recover the Virial theorem obtained in [5,7, 10]: 

±tf = 2Nk B (T-T c ). 
We can define an effective diffusion coefficient 



Deff(T) = 



k B T 



T 



(28) 
0), we 

(29) 
(30) 



so that (r 2 ) = 4D eff (T)t 
grated into 



(r 2 )o- Eq. (29) can be inte- 



I(t) = 1(0) 



ANk t 



■(T-T c )t. 



(31) 



For T > T c , the system evaporates and for T < T c , it 
experiences gravitational collapse in a finite time 



end 



^(0) 



ANk B (T c -T\ 



(32) 



At t = tend, the moment of inertia / = 0, implying that the 
system forms a Dirac peak containing the whole mass M. 
For T = T c , the moment of inertia is conserved (see [5, 7] 
for more details and explicit solutions). 

The above results are very general and do not depend 
whether the system is axisymmetric or not. It is interest- 
ing that we can obtain these qualitative behaviors without 
being required to solve the SP system (l)-(2). 
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Virial theorem for 2D point vortices. We now 

consider a single species system of point vortices in two- 
dimensional hydrodynamics. It has been shown in [18] 
that the Gibbs canonical distribution is reproduced by 
the statistical equilibrium state of a system of Brownian 
point vortices whose iV-body dynamics is defined in terms 
of stochastic equations (instead of deterministic equations 
for usual point vortices [19]). This is the counterpart of 
the system of self-gravitating Brownian particles consid- 
ered previously. In the mean field approximation, the evo- 
lution of the smooth vorticity field u/(r,t) is governed by 
the Fokker-Planck equation [18]: 



^ = V ■ [D (Vw + /3 7 wVV>e//)] 



-Aifj = u, 



(34) 



where 7 is the circulation of a point vortex, /3 = l/(fcsT) 
is the inverse temperature (which can be positive or neg- 
ative [20]) and ip e ff(r,t) = ip(r 7 t) + \Slr 2 is the rela- 
tive stream function in the rotating frame. The Lya- 
punov functional associated to this equation is the Massieu 
function (free energy) J = Sb — /3E e ff where Sb = 
— J — In — dr is the Boltzmann entropy of point vortices 
and E e ff = 5/ utipdr + ^fl J ujr 2 dr is the effective en- 
ergy. Using Eqs. (33)-(34), we easily derive the canonical 
ii-theorem 



J = I — (Vcj + ^wVV'e//) 2 dr > 0. 



At equilibrium, J = 0, implying 
k B T. 



1 



-VOJ + UlVlpeff = 0, 



(35) 



(36) 



which is formally similar to the condition of hydrostatic 
equilibrium (6) for self-gravitating systems with a "pres- 
sure" p = ^^-Lo [11]. Integrating Eq. (36), we obtain the 
mean field Boltzmann distribution 



w(r) = Ae-^°ff ir) . 



(37) 



This distribution is dynamically stable with respect to the 
Fokker-Planck equation (33)-(34) iff it is a maximum of 
Massieu function at fixed circulation. In that case, it rep- 
resents the statistical equilibrium state in the canonical 
ensemble [18]. In the microcanonical ensemble, the sta- 
tistical equilibrium state maximizes the entropy at fixed 
circulation, energy and angular momentum [21,22]. The 
critical points (canceling the first order variations of en- 
tropy) are also given by the mean field Boltzmann dis- 
tribution (37). Therefore, all the results obtained in the 
following at equilibrium, that exclusively use the form of 
Eq. (37) are valid both in the canonical and in the micro- 
canonical ensembles. 



We now derive the equilibrium Virial theorem of point 
vortices. Taking the scalar product of Eq. (36) with r and 
integrating over the whole plane, we obtain 



k B T 

7 



Vw-rdr + V + fi£ = 0, 



(38) 



where we have introduced the angular momentum L = 
J lot 2 dr (similar to the moment of inertia for material 
particles) and defined the "Virial of point vortices" V = 
/ uir ■ Vtp dr by analogy with the Virial of the gravitational 
force. Using the expression of the stream function ip(r) = 
-^/o;(r')ln|r 



find that V = - 



477 



r' I dr' in an infinite domain, we easily 
. Integrating the first term of Eq. (38) 



(33) by parts, we obtain the equilibrium Virial theorem 



r 2 

-2Nk B T - — 
4tt 



m = 0. 



Introducing the mean field critical temperature 

N 7 2 



knTr. = -- 



8- 



it can be rewritten 
1 



-flL = Nk B (T-T c ). 



(39) 



(40) 



(41) 



As explained previously, this relation is valid both in the 
canonical and in the microcanonical ensembles. This re- 
lation was previously derived by Pointin & Lundgrcn [23] 
by a very different method. We show here that it can be 
interpreted as the Virial theorem for a point vortex sys- 
tem. For L = or Q = 0, we deduce that an equilibrium 
state can possibly exist in an infinite domain only at the 
critical temperature T = T c . 

The out-of-equilibrium Virial theorem for Brownian 
point vortices is obtained by taking the time derivative of 
the angular momentum and using the mean field Fokker- 
Planck equation (33)- (34). After straightforward algebra, 
we obtain 



k B T 
4.D7 



1 



L + - flL = Nk B (T - T c ). 



(42) 



Solving this equation, we find that the angular momentum 
evolves like 



L(t) 



L(0) 



2Nki 
O 



(T-T c ) 
2Nk B 



+ - 



(T - T c ). 



Let us introduce the new critical temperature 

flL(0) 



T = T C 



2Nk t 



(43) 



(44) 



This critical temperature depends on the initial value of 
the angular momentum and on the angular velocity. It 
can be written 



To(fi) =T C 1 



n 



(45) 
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with 



2„,2 



4ttL(0) ' 



(46) 



We can define an effective diffusion coefficient 



D ef f(T)=D 1 



T 



(53) 



In terms of the critical temperature (44), the evolution of so that 
the angular momentum is given by 



grated into 



2 ) = AD eff (T)t + (r 2 } . Eq. (52) can be intc- 



L{t) = 



2Nki 



(T - T c ) + (To - T)e 



For small times, we have 



L(t) ~ L(0) 



4L» 7 iV 
T 



(T - T )t, (t 0). 



(47) 



(48) 



In the following, we take 7 > (the case 7 < can be 
treated symmetrically). 

(i) Let us first assume ft > 0. In that case, To > T c . 
If T > 0, then L{t) -> ^^(T - T c ) for t -> +00 (it 
increases for T > T and it decreases for T < T ). The 
system is expected to reach an equilibrium state. If T < 0, 
we need to consider three cases: (a) if T < To, the angular 
momentum Lit) increases and tends to +00 as t — > +00. 
Since L(t) = T(r 2 )(t), the system of Brownian vortices 
"evaporates", (b) If T > To, the angular momentum L{t) 
decreases and becomes zero in a finite time 



tr 



k B T 



In 



T-T c 
T-T 



(49) 



This corresponds to a finite time collapse. At t = t en d, 
the angular momentum L — 0, implying that the system 
forms a Dirac peak containing all the point vortices, (c) If 
T = T , the angular momentum remains constant for all 
times: 



L(t) 



2N 



k B (T -T c ) = L(0). 



(50) 



(ii) Let us now assume < 0. In that case, To < T c , If 
T > 0, the angular momentum L(t) increases and tends 
to +00 as t — > +00, leading to an evaporation process. If 
T < 0, we need to consider three cases: (a) if T < T c , 
then L(t) -> (T - T c ) for i -> +00 (it increases for 
T < T and it decreases for T > T ). The system is 
expected to reach an equilibrium state, (b) If T > T c , the 
angular momentum L(t) decreases and becomes zero in a 
finite time given by Eq. (49). This corresponds to a finite 
time collapse. At t — t ene [, the angular momentum L = 0, 
implying that the system forms a Dirac peak containing 
all the point vortices, (c) If T = T c , then 



L(t) = L{0)e~ 



(51) 



The angular momentum decreases and tends to zero for 
t — > +00, implying that the system forms a Dirac peak 
containing all the point vortices in infinite time. 
For Vi = 0, the Virial theorem (42) reduces to 



4TT , 
L=—{T~T C ). 



(52) 



L(t) = L(0) 



4DNj 
T 



(T - T c )t. 



(54) 



For T > and for T < T c (i.e. > (3 C ), the angular 
momentum increases and tends to +00 for t — > +00 so the 
system evaporates. For T c < T < (i.e. (5 < (3 C ), the 
angular momentum decreases and becomes zero in a finite 
time 



tend — 



-L(0)T 



4T>7V 7 (T - T c 



(55) 



At t = t enc i, the angular momentum L = 0, implying that 
the system forms a Dirac peak containing all the point 
vortices. For T = T c , the angular momentum is conserved. 

Conclusion. — In this Letter, we have shown that the 
Virial theorems for 2D self-gravitating Brownian particles 
and point vortices take a very simple form. While the out- 
of-equilibrium Virial theorems Eqs. (12), (13), (16) and 
(42) are only valid in the canonical ensemble (for Brownian 
systems), the equilibrium Virial theorems Eqs. (10), (17) 
and (41) are valid both in canonical and microcanonical 
ensembles. To make the discussion straightforward, and 
emphasize the basic results, we have considered a single 
species system and made the mean field approximation. 
We shall relax these assumptions in a companion paper. 
Interestingly, we find that the exact Virial theorems keep 
the same form as the mean field Virial theorems except 
that the mean field critical temperatures are replaced by 
the exact critical temperatures [10]: 



JV 



knT r 



k R T r 



—Y y 

AN ^ ^ 



a — l oc'^a 



N 



8irN 



EE 

a — l q'^ql 



la To 



(56) 



(57) 



In the single species case, we get fcsT c = [N — l)Gm 2 /A 
and k B T c = —(N — l)7 2 /(87r) and they only differ from 
the mean field critical temperature by a term of order 
1/N that tends to zero as N — > +00 (recall that the mean 
field approximation becomes exact for N —> +00). The 
mean field Fokker-Planck equations Eqs. (l)-(2) and Eqs. 
(33)-(34) with rotation have never been studied so far and 
the present preliminary results (that have been obtained 
without solving them explicitly) show that they have a 
very rich behavior. This can stimulate further research 
from physicists and applied mathematicians. 
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